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ABSTRACT 



We study the three dimensional fundamental-adjoint SU{2) lattice gauge 
theory at non-zero temperatures by Monte Carlo simulations. On an 8^ x 2 
lattice, at Pa = 1-1, where Pa is the adjoint coupling, we find no evidence 
of any transition at the location of a previously known bulk phase transition 
around /? = 1.33. Moreover, the deconfinement transition at /3a = 1-1 occurs 
at P = 1.20 and is of first order for Pa > 1.1, thus implying a change of 
universality class from that of the Wilson action at Pa = 0. Computations of 
the plaquette susceptibility and the temporal and spatial Polyakov loops on 
8^ X 4 and 16^ x 8 lattices at /3a = 1-1 further support these conclusions and 
suggest that the previously claimed bulk transition around P = 1.33 is, in 
fact, the first order deconfinement transition. Simulations at larger Pa and 
the measurements of the mass gaps from the correlation functions of temporal 
and spatial Polyakov loops also confirm the temperature dependent nature 
of the above transition. The consequences of our results on universality are 
discussed. 



2 



1. INTRODUCTION 



One of the most outstanding problem in physics is to understand the 
mechanism for confinement and finite temperature deconfinement in quan- 
tum chromodynamics (QCD) or SU{N) gauge theories in general. Colour 
confinement implies that the physical states in the Hilbert space are colour 
singlets. Non-abelian gauge theories at zero temperature are believed to have 
this confining property. Regulating a compact gauge theory on the lattice 
with cut-off a~^, where a is the lattice spacing, confinement can be analyt- 
ically established in the strong coupling region. Creutz|l|] showed that no 
phase transition takes place at zero temperature for the usual Wilson ac- 
tion of the SU{2) or SU{3) lattice gauge theory as the continuum limit is 
approached, i.e., these theories have a smooth (3 function for all values of 
the cut-off, with a zero only when the cut-off goes to infinity. In contrast, 
the phase diagrams of lattice 50(3), SU{4) and SU{5) models have bulk 
transitions separating the corresponding strong and weak coupling phases. 
However, Bhanot and Creutz|0] showed that the mere existence of a bulk 
phase transition in lattice theories does not necessarily imply a loss of con- 
finement. As the lattice action of a given continuum theory is not unique, 
one could instead consider an extended lattice theory belonging to the same 
universality class. If the phase diagram of the theory in the extended cou- 
pling space allows a smooth continuation around the bulk singularities then 
the theory will retain its strong coupling confining property as the cut-off is 
removed . They illustrated this by considering an extended Wilson action 
defined by: 



1 - ^TtaUp 



Here F,A denote fundamental and adjoint representations respectively. Up 
is the ordered product of the four directed link variables U^{n){E SU{2)) 
which form an elementary plaquette. The sum over P denotes the sum over 
all independent plaquettes of the lattice. As found in Ref. 0, this model 
has a rich phase structure, which is displayed in Fig. 1. At /3a = 0, it 
reduces to the usual Wilson action. Along the /3 = axis, it describes the 
SO (3) model with a first order phase transition at ~ 2.5. At /?a = cxo, 
it describes the Z2 lattice gauge theory with a first order phase transition 
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at (3^ = iln(l + V2) ~ 0.44. From the p = and Pa = oo axes, the 
above two bulk transitions extend into the Pa) plane, meet at a point 
A at (0.5,2.4) and then continue as another line of bulk first order phase 
transitions. The fact that the latter line ends at a point B at (1.5, 0.9) in the 
phase diagram allows one to bypass all the bulk singularities without losing 
confinement in SU{2) gauge theory. This model therefore retains its strong 
coupling confining property also in the continuum and is expected to be in 
the same universality class as the Wilson action {Pa = 0). 

Comparing the naive classical continuum limit of Eq. (1) with the con- 
tinuum SU (2) Yang Mills action, one obtains 

1/gl = P/4 + 2Pa/3 (2) 

Here Qu is the bare coupling constant of the continuum theory. It is conve- 
nient to characterise the model with Qu and 6, defined by 6' = tan~^/?^//3. 
The two loop /3-function of this theory |^ tells us that theories characterised 
by different ^'s approach the universal critical fixed point g* = in the 
continuum limit, with the flow governed by the asymptotic scaling relation. 



a = l/A{e) exp 

where 



1 



2Po9l 
A(0) 5n^ 6tan9 



9u , (3) 



(4) 

A(^) 11 (3 + 8tan^) ^ ^ 

Here Po and Pi are the first two coefficients of P function for SU{2) gauge 
theory. Eqs. (3) and (4) show that the continuum limit of the theory lies 
at 5'J^=0 and that 6 is an irrelevant coupling which changes only the scale of 
the theory. Up to two loops this illustrates that the extended Wilson action 
described by Eq. (1) is in the same universality class as the Wilson action. 

Usually finite temperature lattice field theory simulations are carried out 
on N^. X A^^ lattices, with iVg- ^ Nj. and the inverse physical temperature 
given by Nra{gu)- Here a{gu) is the lattice spacing at the coupling g^, given 
in the asymptotic scaling region by Eqs. (3-4). The corresponding physi- 
cal volume is V = N^a^. If denotes the transition temperature of the 
deconfinement phase transition and is thus also a physical observable, one 
gets the scaling relation oc in the limit a — 0. Together with Eq. 
(3), this yields the scaling law for the critical coupling g^iNr) for the lattice 
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with its temporal extent A^,-- If N^j is not very large compared to A^t, then 
one expects corrections to this scaling relation which are, however, small: 
The leading effect is still the exponential scaling of A^^ with the correspond- 
ing critical coupling q^INt). Furthermore, the relation between A(^) and 6, 
given in Eq. (4), is obtained in leading order perturbation theory. Correc- 
tions to this relation are scheme dependent and not well known but they 
too must become vanishingly small for sufficiently large lattices. Note that 
as the temporal lattice size becomes infinitely large, i.e., as the continuum 
limit is approached, the critical coupling moves to zero: g^i^r — > c>o) 0. 
Thus in the language of finite size scaling theory one can characterise the 
deconfinement phase transition as one with a critical coupling g* = and 
an effectively logarithmic scaling of the critical coupling with the temporal 
lattice size A^^. 

The finite size scaling behaviour of the shift of the critical coupling (or 
temperature) of a bulk transition [^], on the other hand, is governed by 

[g:{N) - g:{oo)] ^ , (5a) 

and the height of the peak of the corresponding specific heat, which in our 
case is the plaquette susceptibility, should scale according to 

Xmax^N'^' . (56) 

Here 6 is the shift exponent, A^ is the typical size of the system : A^ = 
(N^Nt-Y^^. g^{N) and (^^(00) are the critical couplings for the finite A^'^ 
and the infinite lattices respectively and the location of Xmax can be used to 
define the critical coupling, g'^{N). For a second order bulk phase transition, 
5 = 1/z/ and Us = a/z/, where a and u are the critical exponents of the 
plaquette susceptibility and the correlation length respectively. For a first 
order transition, uj is equal to the dimensionality of the lattice d (=4, in our 
case). The bulk transitions in the phase diagram in Fig. 1 must exhibit 
scaling relations similar to Eq. (5-a,b) for various 6 on the transition line 
AB. As the lattice size is increased, the bulk transition should move towards 
the non-zero g^iO), with a rapid increase in the plaquette susceptibility, and 
remain anchored there as the lattice becomes infinitely large. Thus unlike the 
finite temperature deconfinement transition, the bulk transition is a spurious 
transition not seen by the continuum limit of the theory at ^ 0. 

Recently we argued that the phase diagram in Fig. 1 is incomplete. 
Since the diagram is largely based on Monte Carlo simulations on small 
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finite lattices A'^^ {Ncr = 5-7), it must necessarily have a deconfinement phase 
transition at sufficiently small coupling Qu- Monitoring the order parameter 
(L) for the deconfinement phase transition on A^^ x A^^ lattices, for N„ = 8- 
12 and Nt- = 4 and 6, in order to understand the phase diagram better, we 
found that the known second order deconfinement transition at (3a = 0.0 
moves into the phase diagram as Pa is increased and joins the ffist order 
bulk transition line around Pa — 1-0 as shown by the dotted line QR in 
Fig. 1. This coincidence of the two transitions persists till /5a = 1.5 up to 
which the simulations were performed. Moreover, there was no evidence for 
two separate, bulk and deconfining, transitions. The fluctuations in Polyakov 
loops and plaquettes were correlated and their discontinuities occurred at the 
same coupling. Increasing the temporal size of the lattice resulted in a small 
shift in Pc and the transitions still remained coincident and discontinuous. 
While the critical exponents were in agreement with the Ising model, as 
predicted by universality [|], for Pa < 0.9, there was strong evidence that 
the transition is of ffist order for Pa > l-l- The histograms of |L| showed 
a definite two peak structure with the dip between the two peaks increasing 
with No- a.t Pa = l-l and a co-existence of two phases with a tunnelling time 
(> 10^) at Pa = 1.5. 

The above joining together of the deconfinement transition line with the 
previously claimed bulk transition line leads to a curious paradox because as 
discussed above, both their natures as well as scaling behaviours are different. 
We therefore proposed the following three possible scenarios in our previous 
paper to explain the paradox : 

A] The joining of the deconfinement transition line with the bulk transition 
line is accidental on the lattices used. They will decouple on larger (or 
smaller) lattices. 

B] The entire transition line corresponds to deconfinement phase transi- 
tions only and the previous identification of it as a bulk transition line 
is, in fact, incorrect. 

C] The whole transition line is due to bulk transitions, implying that the 
theory has no confinement in the continuum limit. 

In this paper we try to resolve the above paradox, at least partially, by 
studying the model on 8^ x 2, 8^ x 4 and 16^ x 8 lattices. The spread in A^^ 
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was thus maximised to accentuate the scahng behaviour of the deconfinement 
phase transition. On the 8^ x 2 lattice at (3a = 1-1 we find a deconfining tran- 
sition at /3 = 1.20 with no evidence of any transition around (3 — 1.33, where 
the bulk transition was expected. By also measuring the spatial Polyakov 
line {L'^) in addition to the temporal Polyakov line (L^), we attempted to 
increase the Nj. range further in the direction of asymptotic scaling region. 
The simultaneous measurements of these on the same lattice have a major 
consequence for the bulk phase transition hypothesis; since the volume and 
the geometry of the lattice is unchanged, one expects a bulk transition to 
occur at the same coupling (3 for both of them. We find that the transitions 
in and are located at different couplings for each of the lattices inves- 
tigated. Although we do find the plaquette susceptibility to peak in a range 
of P close to the expected first order bulk transition at Pa — 1-1, the height 
of its peak is found to decrease as the lattice volume increases. Further, the 
deconfinement phase transition at (3a = ^-^ appears to be a first order transi- 
tion, whose strength increases rapidly as one increases (3a'- clear coexistence 
of two-states was seen for Pa = 1-4 and 1.5. The scaling behaviour of the 
transition with Nj- as well as the computations of mass gaps obtained from 
our measurements of the correlation functions of the temporal and spatial 
Polyakov lines again confirm the deconfining and first order nature of the 
above transition. 

The organisation of the paper is as follows: In section 2 we describe the 
observables used to monitor the phase transition. Section 3 is devoted to 
the description of our Monte Carlo data and its analysis. Our results and 

their consequences on universality are discussed in the next section. Finally 
we conclude the paper with a summary and a listing of the still unanswered 
questions. 

2. THE OBSERVABLES 

In order to understand the nature of the phase transitions around the 
transition points of the extended action, defined by Eq. (1), we study the 
expectation values of the temporal and spatial Polyakov loops and the be- 
haviour of their correlators in the two phases. As discussed below, these 
observables and their correlators characterise the phases on the two sides of 
the deconfinement transition. We also studied the expectation value of the 
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plaquette, given by 



(P) = ^-^r>{TrFUr) .g. 



and the corresponding susceptibility, 

XP = QN^NAiP') - (P)') (7) 
Here the expectation value of an observable O is defined by 

(O) = Z-^ In dU^{n)0 exp -S (8) 

with S as the extended action given by Eq. (1) and Z as the corresponding 
partition function. 

Due to periodic boundary conditions in the temporal direction, the SU{2) 
gauge theory has a Z{2) invariance corresponding to the center of the gauge 
group [0. Under this symmetry 

Uo{n,To) ^ zUo{n,To), Vn & fixed tq , (9) 

with z G Z{2). Here Uo{n,T) is the time like link at the lattice site {n,T). 
Under the above transformation, the Polyakov loop, which is the order pa- 
rameter for the deconfinement phase transition and is defined as 

1 1 
L'^jr,T.L'i^) = --Y.Trl[Uo{fi,r) , (10) 

transforms as follows: 

^ z U. (11) 

The expectation value of the Polyakov loop can be shown 0, |^ to be a measure 
of the free energy of an isolated quark. Thus the deconfinement transition 
corresponds to spontaneous breaking of the Z{2) symmetry, separating the 
Z{2) symmetric low temperature (confined) phase from the high temperature 
broken (deconfined) phase. Since due to tunnellings between the two Z[2)- 
vacua the expectation value of this order parameter on a finite lattice is 
zero, we chose, as usual, to look at its absolute value and the corresponding 
susceptibility, defined by 

X = Nl{{Ll)-{\U\f) (12) 
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On an infinite lattice the critical behaviours of the Polyakov loop, its 
susceptibility and the the correlation length of the theory are governed by 
the critical exponents /?, 7 and i/, defined as : 
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(13) 
(14) 
(15) 

One obtains these exponents from Monte Carlo simulations by simply fitting 
the order parameter ^ |10| and/or by using the finite size scaling theory|||]. 
According to the latter, the maximum of the susceptibility, defined in Eq. 
(12), is expected to grow like 

XmaxOCiV,- , (16) 

where a; = 7/z/ for a second order transition and is equal to spatial dimension 
of the system (3 in our case) for a first order transition. 

In addition, we also define averaged spatial Polyakov loops by 

L'' = \j:K , (17) 



with 



E 

n_Lz 



Tr n Uoifi,no) . (18) 



The summation over i is over the three spatial directions and i are the cor- 
responding unit vectors. Due to the periodic boundary conditions in the 
spatial directions also, L'^ can again be thought of as an order parameter for 
the deconfinement transition now for a lattice of spatial volume N^N^- and 
at an inverse temperature of iVg-a. Due to the larger effective "temporal" 
extent, the order parameter L°" can be considered to probe a region closer to 
the asymptotic scaling region compared to that probed by L'^ . Of course, the 
effective 3-volume, N^Nr, seen by L'^ , being smaller than the corresponding 
spatial volume, A*"^, for L'^, finite size effects are relatively more in this case. 
However, as discussed in the next section, on the lattices concerned these 
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corrections only introduce uncertainities of the order of a few per cent in the 
final results discussed in this paper. Therefore, for our purpose of exploring 
the transition region, spatial Polyakov loops are as good an order parameter 
as temporal Polyakov loops. Moreover, due to the averaging over the three 
spatial directions, for Nt- > the statistics for L'^ is also relatively better. 

The most important advantage of measuring both the order parameter L'^ 
and L'^ turns out to be the ability to distinguish between the three scenarios 
mentioned earlier. Let us therefore consider them again in light of these 
observables. 

A] If both bulk and deconfining transitions are present and overlap for 
the lattice sizes investigated in Ref. 0, then the effective increase in 
the temporal extent by studying L'^ may allow us to see the deconfine- 
ment phase transition to shift away. L'^ and L'^ should show critical 
behaviours at /3J and respectively with > At /? < /5J both 
L'^ and L*^ should show confinement. In the region j3l < j3 < j3", 
will be deconfined while L" should still be in the confined phase. For 
(3 > (3^, both L'^ and L'^ should show deconfinement. A presence of a 
bulk phase transition may show up as a non-analyticity in L'^ or in the 
corresponding mass gap obtained from L'^-correlations to be defined 
below. 

B] If there is only a deconfinement transition then the L"^ or the corre- 
sponding mass gap will be smooth at /3^. Moreover, in the asymptotic 
scaling region the difference — should decrease and approach a 
constant as we increase the lattice size holding Nu/Nj. constant while it 
should increase if the lattice size is increased only in the A^o--directions. 

C] If there is only a bulk transition present then /^^T ~ : both L'^ and 
L°" should become significantly different from zero at about the same 
beta, especially since they are evaluated on a lattice with identical 
4-volume and geometry. Also the average plaquette should show a 
critical behaviour roughly at the same coupling. Unlike scenario [B], 
the differences in the critical couplings should vanish in both the limits 
discussed there. 

Another useful observable, mentioned already above, is the mass gap 
obtained from the correlations of the Polyakov loops. For this purpose, we 
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define sums of temporal and spatial Polyakov lines over the corresponding 
orthogonal planes: 

^>.) = ^E^^(^) (19) 

^ n±j 

^rK) = T7V E Lnno,n) (20) 

-^^a-^^T riQ, 

n_Li ,n_Lj 

These can be thought of as the zero momentum projections of the Polyakov 
lines in 2+1 dimensions. In our simulations, we measure the zero momentum 
correlation functions of both the temporal and spatial Polyakov loops given 
by: 

T--(--\\n\) = {r^''\0)r^''\n)) , (21) 

Here the indices on have been suppressed. Using periodic boundary 

conditions in the temporal and spatial directions and introducing the eigen- 
states of the transfer matrix as the intermediate set of states one obtains 



\n\ 



r^^(-'^)(|n|) = Z-i5:(/|/:^('^)(0)|m) exp[-/i;,('^)|n|]exp[-/x[(^)(iV.- 

(22) 

In the large distance limit, ignoring the higher excited states of the transfer 
matrix and summing over only the lowest states, one can show that 

r--(--)(|^|)=^-H2|exp-/x"("Vl+exp-;u"(")(Ar, -|n|]} , (23) 
with = Z^^(0|£^*-'^''|1)^. /i^^'^-* are the mass gaps which can be shown[|rT], 



1% , p!3| , to correspond to the physical mass, namely, string tension, and 
the tunnelling mass in the confined and deconfined phases respectively. This 
is because of the expectation of a linearly rising colour averaged potential 
between a heavy quark and antiquark in the confined phase. In the de- 
confined phase, the potential may become of the Debye-screened Coulomb 
form. However, on a finite lattice, the tunnelling between the two degener- 
ate Z(2)-vacua, | 0+) and |0_) yields a symmetric |0s) ground state and an 
antisymmetric first excited state \0a), leading to a small "tunnelling mass 
gap". In the infinite volume limit, the tunnelling mass goes to zero and the 
symmetry is spontaneously broken. Thus this mass gap also acts as an ef- 
fective order parameter for the spontaneous breaking of the Z(2)-symmetry 
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and deconfinement. Its main advantage over the Polyakov loop is its direct 
relation with the string tension in the confined phase. Thus its becoming 
vanishingly small across a phase transition is a clearer indication of it be- 
ing a deconfinement phase transition. Of course, the string tension will also 
decrease across a bulk phase transition but its vanishing on the other side, 
especially as the lattice size grows to infinity, is very unlikely. 

3. DATA AND ANALYSIS 

Our Monte Carlo simulations were done mostly at (3a = 1-1 on 8^ x 2, 
8^ X 4 and 16^ x 8 lattices, using the Metropohs et al. algorithm. On the 
8^ X 2 lattice we also performed simulations at larger Pa to confirm the order 
of the phase transition and to estimate the shift in critical f3. Additional 
simulations were also made on 10^ x 2 and 12^ x 2 lattices to look for the 
finite size scaling behaviour and the critical exponent uo, defined in section 2. 
Typically 1-2 xlO^ iterations were done on all lattices with about a factor 
two less statistics on the biggest lattice. For the 8^ x 4 and 16^ x 8 lattices, 
we also measured the spatial Polyakov loops. Measurements of temporal and 
spatial zero momentum correlators were done on the 16^ x 8 lattice both in the 
confined and deconfined regions. The distribution function of the Polyakov 
loop is a good tool for identifying the critical points to a reasonably good 
accuracy, especially for a first order transition on larger lattices. In the 
confined phase one expects this function to be a gaussian peaked around 
zero, while in the deconfined phase of the SU(2) gauge theory one expects 
two symmetric gaussian peaks located away from zero. At the transition 
point of a first order phase transition, one would have an additional peak 
of the same height located at zero. One can alternatively look at the order 
parameter, \L\, and correspondingly obtain a two peak structure. Even when 
the transition is of second order, or the precise location of the transition point 
is not found, these distribution functions can be used to set the lower and 
upper limits on (3(.- For our semi-quantitative purposes of distinguishing the 
shifts due to the two different scaling behaviours of the deconfinement and the 
bulk phase transition on 8^ x 4 and 16^ x 8 lattices, it was usually enough to 
study the behaviour of the histograms of the temporal and spatial Polyakov 
loops at various /3-values. 
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3.1 = 2 



Variation of the temporal size of the lattice should shift the deconfinement 
phase transition much more than the bulk phase transition, allowing one to 
test the scenario [A] of the previous section. Already in Ref. ||^, a change 
from A^^ = 4 to 6 was attempted but the results were inconclusive as the shift 
in the critical coupling was too small and only one transition was observed. 
From the scaling law for the deconfinement phase transition, one expects a 
much larger shift by going to smaller N-r, although it usually takes one in the 
strong coupling region. Since our aim here is only to see whether there are 
two separate transitions, it is perhaps not so important that one does not 
work in the scaling region; hopefully, such qualitative aspects are not altered 
by the continuum limit. 

Figs. 2a and 2c show the average temporal Polyakov line (L^) and the 
average plaquette as a function of the coupling /5 for = 1.1 on an 8^ x 2 
lattice. The later value was chosen because our earlier work0 suggests a 
weak first order phase transition at {(3, (3a) = (1.33, 1.1) on an 8'^ x 4 lattice 
which could be coincident with the bulk transition seen in Ref. [0 . The 
corresponding susceptibilities are displayed in Figs. 2b and 2d. Both the 
Polyakov line and the plaquette have a sharp and large jump around /? = 
1.20 and correspondingly the respective susceptibilities exhibit a sharp peak 
there. The above behaviour of the Polyakov line and its susceptibility clearly 
indicates the presence of a deconfinement transition at f3c = 1-20. Note, 
however, that the behaviour of the Polyakov line and the plaquette is smooth 
in the region around (3 = 1.33 where the bulk transition was expected. For a 
better comparison a magnified view of the Polyakov line and the plaquette in 
the regions around the deconfinement phase transition and around (3 = 1.33 
is shown in the Figs. 3a - 3d. In order to show a bigger region around 
f3 = 1.33, ranges of both x and y axis have been enlarged by a factor of four 
in Figs. 3b and 3d as compared to those of Figs. 3a and 3c. The remarkable 
fiatness of the data in Figs. 3b and 3d, together with the absence of any 
structure in the susceptibilities in Figs. 2b and 2d, make it clear that there 
is no trace of any transition in the region around jS = 1.33. In view of the 
different scaling laws for the deconfinement and bulk phase transitions and 
the additional results from Ref. , indicating a very small shift in f3c in going 
from an 8^ X 4 lattice to a 12^ x 6 lattice, this relatively large shift makes it 
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difficult to sustain the accidental coincidence hypothesis in scenario [A]. Of 
course, the lattices used may not be large enough for various corrections to 
the scaling behaviour to be small. Nevertheless, it would be remarkable, if 
not impossible, if they are similar in magnitude and sign for the two unrelated 
scaling forms. 

The large shift in the critical coupling in going to A^^ = 2 from A^^ = 4 
and 6, and the fact that the order parameter (|L|) shows a spontaneous 
breakdown of the Z{2) symmetry at the transition point are indicative of the 
deconfining nature of the phase transition. Moreover, the histograms of L, 
on 8^ X 2, 10^ X 2 and 12^ x 2 show a three peak structure at the transition 
point which, however, is not very pronounced at (3a = l-l- Correspondingly, 
the determination of the critical exponent u does not fix the order of the 
phase transition uniquely. This is similar to the A^^ = 4 results in Ref. [^]. 
It was found there that the above transition was a relatively weak ffist order 
transition at /3a = 1.1 but for (3 a > 1-1 the discontinuity increased with an 
increasing Pa- Having thus verified that the order of the phase transition 
does not change as A^^ is changed to 2 from 4, we simulated the theory at 
Pa = 1-4 and 1.5 on the above lattices with A^^ = 2 to see if its strength 
increases. It does indeed increase and we were even able to conffim the ffist 
order nature by determining the critical exponent u. The evolution curves 
for both the plaquette and the Polyakov loop in the respective critical regions 
at Pa = 1-4 are plotted in the Figs. 4a-4c and the corresponding histograms 
for the latter are shown in Fig. 5. One sees that they all suggest a clear 
ffist order deconfinement phase transition. One does observe a sharpening of 
the transition as the volume is increased. It may be interesting to note that 
the average plaquette also shows a discontinuity and fluctuations which are 
totally correlated to that of L, which is again in agreement with the similar 
observation on A^^ = 4 lattice. Individual runs around the critical coupling 
were used to locate the maximum of the Polyakov susceptibility. This maxi- 
mum value of X was in turn used to compute the exponent uj, defined earlier 
by Eq. (16). On the 12^ x 2 lattice, a large tunnelling time of (^30,000-40,000 
sweeps) was observed at Pa = 1-4, and hence a relatively higher statistics 
run of 3 X 10^ was made at the transition point to ensure a better sampling 
of both the phases. The critical couplings and the estimated values of the 
exponent u are given in the Table 1. They are in excellent agreement with 
the finite size scaling prediction for a ffist order deconfinement phase tran- 
sition. Finally, also the simulations at Pa = 1-5 support the above picture 
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and clear three-peak histograms are observed. The critical coupling for the 
8^ X 2 lattice for (3a =1-5 is /? = 0.983 ± 0.001. Comparing this value with 
the corresponding A''^ = 4 result, one sees that the shifts due to a change of 
A^^ from 4 to 2 become smaller as (3a grows larger, but even at (3a = 1-5 the 
shift is nevertheless significantly large enough to be sure that it is non-zero. 
The discontinuity in |L| also appears to increase slightly on the smaller N^- 
lattice when compared at the same (3a- 

3.2 = 4 and 8 

The above results on the Nr = 2 lattices showed clearly the absence of 
any transition at the location of a previously claimed bulk phase transition 
and further showed a first order phase transition at a much smaller (3. The 
latter can be consistently interpreted as a deconfinement phase transition. 
While they make the coincidence of a bulk transition unlikely, one would 
like to establish it as concretely as possible. With such a motivation, we 
simulated the extended action at (3a = 1.1 on 8^ x 4 and 16'^ x 8 lattices. 
Now, we also monitored the L"^ since the corresponding effective 3-volume is 
substantial on these relatively bigger lattices. 

The histograms for and on the 8^ x 4 lattice at four values of 
the coupling (3 = 1.31, 1.32635, 1.4 and 1.65 are shown in Figs. 6a-6d. 
These couplings have been chosen such that the two intermediate couplings 
are close to the corresponding (3(( and (3^. Due to the weakness of the first 
order transition, it was not possible to get a clear three-peak structure of 
the histogram. However, the reasonably flat behaviours of the Lr and L^- 
in the Figs. 6b and 6c show the onset of the deconfining transition and 
suggest that it is first order. Locating the critical point by the criterion of 
maximum susceptibility, /5J was found in Ref. ||^ to be 1.327 which is good 
agreement with the (3^ = 1.32635 of Fig. 6c. Figs. 6a-6d, and our additional 
simulations (not shown in the these figures) just below and above the critical 
couplings, show that these distributions follow the expected pattern of the 
deconfinement transition very neatly: 

a] For (3 < (31, both L^- and L„ are confined, 

b] at /3 = = 1.32635, is critical and L„ is confined, 
c\ 3X (3 = (3^ = 1.4, Lo- is critical and Lr is deconfined, 

d] for (3 > P^, both Lr and Lo- are deconfined. 
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A similarly chosen pattern of histograms on the 16^ x 8 lattice at five 
values of the coupling /?=1.33, 1.35, 1.3508, 1.47 and 1.57 is shown in Figs. 
7a-7d. Again, we observe a confinement-deconfinement pattern in and 
L„ like the one discussed above. Due to the narrower transition region, we 
were unable to get better histograms of Lj. and fix the critical coupling (3^ 
more accurately. However, as is clear from the Fig. 6b, the critical coupling 
/3^(A^T- = 8) is between 1.35 and 1.3508. We have also observed a clear 
signal for the coexistence of two phases at (3 = 1.35, with a discontinuity 
in the Polyakov line which is roughly 0.05. The corresponding value 
at the critical point on the 8^ x 4 lattice was found to be 0.25. Due to 
the fiuctuations in plaquette, its discontinuity could not be estimated but it 
appeared to be smaller than the one observed on the 8^ x 4 lattice as well. 
Such a scaling of discontinuities is expected for a deconfinement transition 
but not a bulk transition. Fig. 8a illustrates the scaling of the Polyakov line 
and the deconfinement temperature at = 1.1 as we increase the temporal 
lattice size and thus approach towards the continuum limit. One can clearly 
see the transition moving towards the continuum critical fixed point 5'J^=0 
with a corresponding reduction in the discontinuity of L^. On the other hand 
the scaling behaviour of the plaquette susceptibility at /?a = 1.1 on the same 
lattices, namely on 8^ X 2, 8^ X 4 and 16^ X 8 lattices, shown in Fig. 8b, 
does not show any bulk phase transition. It should be noted that on the 
largest lattice the peak of the susceptibility is indeed at a bit lower value of 
/3 compared to the location of the deconfinement transition hinted in Figs. 
7-b and 8-a. Such a behaviour would be consistent with the scenario [A] 
where one expects the deconfinement phase transition to move away faster. 
However, the height of the peak of the susceptibility for the 16^ x 8 lattice is 
seen to decrease in Fig. 8b although the bulk volume increases by a factor of 
16. Combined with the expected finite size scaling behaviour of a bulk first 
order in Eq. (5b), this result clearly rules out the presence of bulk transition 
convincingly in this region too. 

Comparing the critical value (3^ = 1.4 for L'^ = -^(8) on the 8'^ x 4 with 
the more precise value Pi = 1.3504 ± 0.0004 for = L(8) on the 16^ x 8 
lattice, one finds that the former is slightly overestimated . If one uses the 
maximum of the susceptibility xdLgl) to define the transition point on the 
8^ X 4 lattice on the other hand, then one obtains (3^=1.348, where the 
Ferrenberg-Swendsen ||15|| histogramming method was used on the configura- 
tions generated at (3 =1.4 to locate the peak. This relatively large difference 
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in the determination of (3c by two different methods for the 8^ x 4 lattice is 
perhaps indicative of the smallness of the effective 3- volume for Lg. Even 
on this small lattice thus the errors in determination of /3c by employing L^. 
as the order parameter are ~ 5 % and a quantitative estimate of (3c on an 
effectively larger temporal lattice is thus feasible. Since we are anyway only 
interested in seeing whether the transition indeed marches towards g* = 0, 
such small shifts are perhaps tolerable to us, considering that it will allow 
us to probe Nr = 16 by studying on 16^ x 8 lattice. With this caveat in 
mind, we see from Fig. 7c that (3c — 1.47 for A^^ = 16, as the Lie-distribution 
can be seen to become very flat at L = on the 16^ x 8 lattice at this cou- 
pling. Furthermore, using the configurations at (3 = 1.47 to locate the peak 
of the |L|-susceptibility, one obtains = 1.468. Considering the increased 
effective 3-volume in this case, such an excellent agreement in these two de- 
terminations is in accord with the expectations. Allowing for a few per cent 
overestimation nevertheless due to the unusual geometry, one still sees that 
the transition is indeed slowly moving towards Qu = 0, although compared 
to the Pa = 0.0 results of Refs. |T^, ^ these results are still far away from 
exhibiting a consistency with the two-loop scaling equations in Eq. (3-4). 

The systematic behaviour of the histograms in both Figs. 6 and 7, sug- 
gesting that L"^ and L"^ start being non-zero at different locations with the 
former doing so always earlier, has important consequences for the scenarios 
[A] and [C] of sect. 2. Since the 4-volume for each figure is held constant, 
and so is the geometry of the lattice, merely one bulk phase transition would 
be inadequate to explain this behaviour of the Monte Carlo data. Thus, 
the scenario [C] is untenable. The strong correlation of the critical coupling 
with only the linear extent of the lattice along which the order parameter 
is defined, taken together with the decrease in the plaquette susceptibility 
at its peak in Fig. 8b and the results of previous subsection for Nr = 2 
lattice, is also inconsistent with scenario [A], thus leaving the scenario B of a 
pure deconfinement phase transition as the only plausible one. Furthermore, 
the above mentioned decrease in the size of discontinuities in the lattice ob- 
servables, namely the Polyakov line and the plaquette, again confirms the 
confining-deconfining nature of the transition. For a bulk transition these 
discontinuities would be expected to remain unaltered or increase, as the 
total volume of the lattice increases. 



17 



3.3 Polyakov Loop Correlations 



As a final check of the nature of the transition along the line QR, which 
shifts strongly with the "temporal" extent of the lattice, as we saw above, 
we have investigated the behaviour of the Polyakov loop correlations and 
the corresponding mass gaps on the 16^ x 8 lattice. Prom the discussion 
above and in sect. 2, one expects to see the mass gaps from the temporal 
and spatial Polyakov loop correlations to decrease strongly at the respective 
critical couplings. Moreover, they should roughly be the same, apart from 
finite size effects, in the phases where both the order parameters show the 
same behaviour. 

The measurements of the zero momentum projected temporal and spa- 
tial correlation functions both above and below the transitions were done 
by recording their values after every 30 sweeps. The individual averages of 
temporal and spatial correlations were taken by summing r^'^(nj) over the 
three directions (j) and rf'^(nj) over (i,j) for fixed value of Uj. The discrete 
rotational covariance was ensured by plotting them for each of the three spa- 
tial directions (i=l,2,3) individually and ensuring that their differences lie 
within error bars. The above averaged values of the temporal and spatial 
correlators are plotted in the Figs. 9a-9b. A gradual flattening of the cor- 
relation functions is clearly seen in both Figs. 9a and 9b, as one increases 
the coupling 13 through the respective 13^- Below the transition they all fall 
rather steeply while a very flat behaviour is evident for very large j3. The 
mass gaps in the conflned and deconflned phases were obtained by fltting 
the correlation function with a single hyperbolic cosine function. Their val- 
ues and the corresponding errors are given in Table 2 and are plotted in 
Fig. 10. The relatively large fluctuations suggest that much higher statis- 
tics may be needed for a precise and reliable determination of these mass 
gaps. Nevertheless, the small value of the mass gap above the respective 
phase transitions, indicated by the corresponding Polyakov line, is a clear 
signal of the Z{2) symmetry being broken spontaneously in the infinite vol- 
ume limit, leading to deconfinement of quarks. The relatively larger finite 
size effects for the spatial correlations are evident in the slower fall-off of the 
corresponding mass gaps, especially in the deconfinement region. The results 
are again in agreement with the conclusions of the previous subsections: one 
has deconfinement phase transitions at the respective (5c. Although the data 
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are somewhat noisy and thus not very conclusive, there is again no evidence 
of any extra first order bulk phase transition. It would have been nicer to 
have the spatial mass gaps at still lower (5 to confirm a continuity in them. 
However, we found a very rapid deterioration of the signal in the correlation 
functions in that region and large fluctuations in the correlation functions 
prevented us from arriving at any values. 
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4. DISCUSSION AND CONCLUSIONS 



Our results in the previous section strongly suggest that the previously 
identified line of first order bulk phase transitions, shown in Fig. 1 by the 
line QR, is in fact a line of first order deconfinement phase transitions. No 
bulk phase transition seems to be present in that region. Since this bulk 
transition and its "shadow" on the /3a = axis have been speculated to have 
an impact on the onset of confinement in the SU{2) gauge theory at zero 
temperature, it is clearly necessary to investigate the phase diagram of Fig. 
1 afresh on much bigger symmetric lattices than used in the past. Similar 
structure is also known to exist for other SU{N) gauge theories, in particular 
for the SU{3) theory as well. These may need to be investigated further to 
investigate their differences. 

A first order deconfinement phase transition for [3^ > 1-1 also implies a 
change of the universality class for the extended action at and above Pa = 
1.1. It has been argued that by integrating out the irrelevant degrees of 
freedom of the SU{2) theory at finite temperature, one can obtain an effective 
theory of the order parameter L. If this theory and the Ising model in 3- 
dimensions, which too has a global Z{2) symmetry for its order parameter, 
have only one fixed point in the space of couplings, then the two are in 
the same universality class. This leads to a prediction of various critical 
exponents for the former. For /3a = 0, these universality predictions were 
explicitly verified by high precision Monte Carlo simulations of Engels et 
al.H. They found f = 0.545 ± 0.030, uj = ^ = 1.93 ± 0.03, u = 0.65 ± 0.04. 
The corresponding values for the Ising model are ^ = 0.516 ± 0.005, ^ = 
1.965 ± 0.005, u = 0.63 ± 0.003. The change of the order of the transition 
(e.g at Pa = 1-4, ~ 3) as a function of an apparently irrelevant coupling Pa 
shows a non-universal behaviour of the deconfining transition, which needs 
to be understood in the renormalisation group picture. It does cast a doubt 
on our understanding of the order of the deconfinement phase transitions 
from those in the corresponding spin models. Of course, one cannot rule 
out the possibility that this change of universality class is a purely strong 
coupling phenomenon; simulations on sufficiently larger lattices will show 
a second order deconfinement phase transition even at the large Pa used 
in this work. It would be very interesting, and illuminating, to check this 
since no such instance of a change of the order of the phase transition in 



20 



going over to the scaling regime is known so far. If confirmed, it would 
underline the importance of checking even such qualitative aspects of full 
QCD thermodynamics in the scaling or the asymptotic scaling regime. 

A non-universal behaviour of the string tension depending on the lattice 
action used has been a question of debate in the past. In the SU(2) lattice 
gauge theory the A-ratios of the Manton and the Heat kernel action with that 
of the Wilson action were found to differ from their theoretically predicted 
values. However, these discrepancies were shown [|^] to be due to higher 
order corrections to the theoretical values which were obtained in leading 
order perturbation theory. By taking the ratios of the physical observables 
in Monte Carlo simulations or by considering the effects of higher order terms 
the above actions were shown to be in the same universality class. In the case 



of the extended SU(2) action of Eq. (1), Bhanot and Dashen |jT9[ found that 
at [3a = 0.9, the A-ratio, I^Extended/ ^wusom obtained from the computations 
of string tension in Monte Carlo simulations, was roughly a factor of four 
higher than its perturbation theoretical value in leading order. Again, this 
was shown p!8|, pO| to be due to the strong coupling effects. In fact, the 



string tension data compared well with its strong coupling expansion pO 
and correspondingly the 3-loop contribution to A ratio was quite large |18 . 

What we find in the change of the order of the deconfinement phase tran- 
sition, however, is a loss of universality in qualitative features, and therefore 
should be considered much more seriously. Moreover, the line of deconfine- 
ment transition we investigated on = 4 is remarkably close to the line 
of constant string tension, obtained from Monte Carlo simulation [1^, at 
(TO? = 0.14. This line lies in the cross-over regionpOf. Thus our results 
and conclusions could possibly be free from any strong coupling artifacts. In 
other words, one need not expect that by going to larger temporal lattices at 
/Sa > 1-0 the scaling exponent u) will drop from 3 to its second order value 
1.97, corresponding to the 3-dimensional Ising model. Taking seriously the 
above coincidence of the lines of the deconfinement transition and the con- 
stant string tension, along with the fact that up to 2-loops of the P function 
physics is independent of the value of 9, the above change of universality 
class indicates a possibility of a new fixed point and some non-perturbative 
effects. However, much more work is necessary to convincingly establish or 
rule out the above interesting possibilities. 

In the case of pure Wilson action, Polonyi and Szlachanyi showed the 
second order nature of the deconfining transition by computing the effective 



21 



action in terms of the order parameter L in the strong coupling limit. It will 
be interesting to compute the effect of adjoint coupling (3a and see how it 
changes the order of the transition. One expects a term with its coefficient 
becoming negative as the adjoint coupling increases above 13a = l-l- 

Let us summarise by stressing that the universality hypothesis in lattice 
gauge theories, in particular of extended Wilson action, should be be looked 
at more critically. It will be interesting to study similar phenomenon with 
the corresponding extended SU(3) Wilson action also. 
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FIGURE CAPTIONS 



Fig.l The phase diagram of the extended SU(2) lattice gauge theory. The 
sohd hnes are from simulations done on a 5^ lattice by Bhanot and Creutz[^ 



The broken lines QR is the finite temperature deconfinement phase transition 
line on the 8^ x 4 lattice. The deconfinement line ST corresponds to 8^ x 2 
lattice. 

Fig. 2 (a)Average temporal Polyakov loop, (b) its susceptibility, (c) average 
plaquette, and (d) its susceptibility vs. /? on an 8^ x 2 lattice at Pa = l-l- 

Fig. 3 Average temporal Polyakov loop and average plaquette as a function 
of on an 8^ X 2 lattice in the regions of the deconfinement phase transition 
(a, c) and the expected bulk transition (b, d). 

Fig.4 Evolution of and (P+0.35) Pa = 1-4 and on (a) 8^ x 2 {p^ = 1.03954) 
(b) 10=^ X 2 = 1.0398) and (c) 12^ x 2 {p^ = 1.04007) 

Fig. 5 Normalised probability distributions of Lr on an 8^ x 2, 10^ x 2 and 
12^ X 2 lattices at Pa = 1.4 and /?c=l. 03954, 1.0398 and 1.04007 respectively. 



Fig. 6 Normalised probability distributions of L"^ and L'^ on 8^ x 4 lattice at 
Pa = 1.1 and (a) P = 1.31, (b) P = P^ = 1.32635, (c) P = P^ ^ 1.4 and (d) 
P = 1.65 

Fig. 7 Same as Fig. 6 but for a 16^ x 8 lattice and at (a) P = 1.33, (b) 
p = p^ = 1.35, 1.3508, (c) P = P^ ^ 1.47 and (d) p = 1.57 

Fig. 8 The nature of the first order transition at = 1.1: (a) The scaling 
behaviour of on 8^ x 2, 8^ x 4 and 16^ x 8 lattices in accordance with 
the deconfinement transition and (b) the scaling behaviour of the plaquette 
susceptibility, ruling out any bulk transition in this region. 

Fig. 9 The (a) temporal and (b) spatial correlation functions on 16^ x 8 
lattices at /?a = 1.1 as a function of distance for various couplings P through 
the respective Pc- 

Fig. 10 The temporal and spatial mass gaps fi'^ and on 16'^ x 8 lattice at 
Pa = 1.1. 
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Table 1 



The values of the critical couplings {PciPa = 1-4) on A^V=2 lattices, the 
finite size scaling exponent lo. The expected value for u; is 3.0(1.97) 
if the deconfining phase transition is first order (second order) . 









8 


1.039540 




10 


1.03980 


3.246(243) [8=^ x 2 : 10^ x 2] 


12 


1.040070 


3.204(184) [12^ X 2 : 8^ X 2] 
3.154(406) [10^ X 2 : 12^ X 2] 



Table 2 



The values of the temporal and spatial mass gaps, /jl^ and //^, at (3a — 1-1 
on 16^ X 8 lattice. 









1.34 


0.33526(01368) 
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0.16618(00568) 
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